DAMAGE SEGREGATION AT FISSIONING MAY 
INCREASE GROWTH RATES: A SUPERPROCESS MODEL 
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Abstract. A fissioning organism may purge unrepairable damage by 
bequeathing it preferentially to one of its daughters. Using the mathe- 
matical formalism of superprocesses, we propose a flexible class of ana- 
lytically tractable models that allow quite general effects of damage on 
death rates and splitting rates and similarly general damage segregation 
mechanisms. We show that, in a suitable regime, the effects of random- 
ness in damage segregation at fissioning are indistinguishable from those 
of randomness in the mechanism of damage accumulation during the or- 
ganism's lifetime. Moreover, the optimal population growth is achieved 
for a particular finite, non-zero level of combined randomness from these 
two sources. In particular, when damage accumulates deterministically, 
optimal population growth is achieved by a moderately unequal division 
of damage between the daughters, while too little or too much division 
is sub-optimal. Connections are drawn both to recent experimental re- 
sults on inheritance of damage in unicellular organisms, and to theories 
of aging and resource division between siblings. 



1. Introduction 

One of the great challenges in biology is to understand the forces shaping 
age-related functional decline, termed senescence. Much current thinking 
on senescence (cf. [68]) interprets the aging process as an accumulation of 
organismal damage. The available damage repair mechanisms fall short, it 
is often argued, because of limitations imposed by natural selection, which 
may favor early reproduction, even at the cost of later decrepitude. One 
line of research aims to clarify these trade-offs by examining the non-aging 
exceptions that test the senescence rule. Of late, it has even been argued that 
negligible [32] or even negative [75] senescence may not be as theoretically 
implausible as some had supposed, and that it might not even be terribly 
rare [5U] . 

Fissioning unicellular organisms have been generally viewed as a large 
class of exceptions to the senescence rule. Indeed, their immortality has been 
considered almost tautological by the principle enunciated by P. Medawar 
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[54j . that individual birth is a fundamental prerequisite for aging. This 
principle has been sharpened by L. Partridge and N. Barton [59], who remark 

The critical requirement for the evolution of ageing is that 
there be a distinction between a parent individual and the 
smaller offspring for which it provides. If the organism breeds 
by dividing equally into identical offspring, then the distinc- 
tion between parent and offspring disappears, the intensity of 
selection on survival and reproduction will remain constant 
and individual ageing is not expected to evolve. 

Recent experiments [H[5lJ[3j[70] have focused attention on the elusive quality 
of the "distinction" between parent and offspring. If aging is the accumula- 
tion of unrepaired damage, then "age" may go up or down. The metazoan 
reproduction that results in one or more young (pristine) offspring, entirely 
distinct from the old (damaged) parent, is an extreme form of rejuvena- 
tion. This may be seen as one end of a continuum of damage segregation 
mechanisms that include the biased retention of carbonylated proteins in the 
mother cell of budding yeast [1] and perhaps the use of aging poles inherited 
by one of the pair of Escherichia coli daughter cells as, in the words of C. 
Stephens [69j , cellular "garbage dumps" . Even where there is no conspicuous 
morphological distinction between a mother and offspring, the individuals 
present at the end of a bout of reproduction may not be identical in age, 
when age is measured in accumulated damage. Whereas traditional the- 
ory has focused on the extreme case of an aging parent producing pristine 
offspring, it now becomes necessary to grapple with the natural-selection 
implications of strategies along the continuum of damage-sharing between 
the products of reproduction. 

Our approach is a mathematical model of damage-accumulation during 
a cell's lifetime and damage-segregation at reproduction that quantifies (in 
an idealized context) the costs and benefits of unequal damage allocation 
to the daughter cells in a fissioning organism. The benefits arise from what 
G. Bell [8] has termed "exogenous repair": elimination of damage through 
lower reproductive success of individuals with higher damage levels. 

For a conceptually simple class of models of population growth that flex- 
ibly incorporate quite general structures of damage accumulation, repair 
and segregation, we analytically derive the conditions under which increas- 
ing inequality in damage inheritance will boost the long-term population 
growth rate. In particular, for organisms whose lifetime damage accumula- 
tion rate is deterministic and positive, some non-zero inequality will always 
be preferred. While most immediately relevant for unicellular organisms, 
this principle and our model may have implications more generally for the- 
ories of intergenerational effects, such as transfers of resources and status. 

One consequence of exogenous repair may seem surprising: If inherited 
damage significantly determines the population growth rate, and if damage 
is split unevenly among the offspring, there may be a positive benefit to 



DAMAGE SEGREGATION 



3 



accelerating the turnover of generations. In simple branching population 
growth models, the stable population growth rate is determined solely by 
the net birth rate. In the model with damage, increasing birth and death 
rates equally may actually boost the population growth rate. This may be 
seen as the fission analogue of Hamilton's principle [10] linking the likelihood 
of survival to a given age with future mortality-rate increases and vitality 
decline, and placing a selective premium on early reproduction. Of course, 
if the variance in damage accumulation is above the optimum, then this 
principle implies that reducing the inequality in inheritance, or decreasing 
birth and death rates equally, would be favored by natural selection. 

2. Background 

Popular reliability models of aging (such as [501 E51 ED] and additional 
references in section 3 of [68] ) tend to ignore repair, while the class of growth- 
reproduction-repair models (such as [321 El CEZ1 E21 [11] and further references 
in section 2 of [68] ) tend to ignore the fundamental non-energetic constraints 
on repair. A living system will inevitably accumulate damage. Damage- 
repair mechanisms are available, but these can only slow the process, not 
prevent it. The repair mechanism of last resort is selective death. Any 
individual line will almost certainly go extinct, but the logic of exponential 
profusion means that there may still some lines survivingQ Following this 
logic one step further brings us to our central question: What benefit, if 
any, would accrue to a line of organisms that could not purge damage, but 
could selectively segregate it into one of its children? 

This mechanism of selective segregation of damage was described by G. 
Bell [8], who called it "exogenous repair". Bell's theoretical insight derived 
from his analysis of a class of experiments that were popular a century 
ago, but have since largely been forgotten. The general protocol allowed 
protozoans to grow for a time, until a small sample was plated to fresh 
medium, and this was repeated for many rounds. Effectively, this created an 
artificial selection of a few surviving lines, selected not for maximum fitness 
(as would be the case if the population were allowed to grow unmolested), 
but at random. In most settings, in the absence of sexual recombination, 
population senescence — slackening and eventual cessation of growth - 
was the rule. While Bell attributed this decline to Miiller's ratchet and the 
accumulation of genetic damage, the same principle could be applied equally 
to irreparable somatic damage. 

One would like to follow the growth of individual cells, but this was in- 
conceivable with the technology of the time. Only lately has such an un- 
dertaking become not only conceivable, but practicable. Recent work [70] 
has followed individual E. coli over many generations, following the fates 

*While this argument is simplest for haploid fissioning organisms, the same logic could 
be applied to the germ line of higher organisms, as in [5]. Sexual recombination, in this 
picture, only facilitates the rapid diffusion of high-quality genetic material, and the bodies 
are rebuilt from good raw materials in each generation. 
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of the "old pole" cell (the one that has inherited an end that has not been 
regenerated) and the "new pole" cell. Whereas individual cells do not have 
a clearly defined age, it makes more sense to ascribe an age to a pole. Un- 
published work from the same laboratory (described in October 2004 at a 
workshop held at the Max Planck Institute for Demographic Research in 
Rostock), as well as the work of [H [51] on Saccharomyces cerevisiae goes 
even further, tracking the movement of damaged proteins or mitochondria 
through the generations, and showing that the growth of the population is 
maintained by a subpopulation of relatively pristine individuals. The demog- 
raphy of damage accumulation in the fissioning yeast Schizosaccharomyces 
pombe has recently been described in [55]. These striking experiments have 
helped focus attention on the elusive nature of the asymmetries that underly 
microbial aging. 

Suppose the "costs" of accumulating damage are inevitable, and that only 
the equality of apportionment between the two daughters may be under evo- 
lutionary control. Of course, perfectly equal splitting of damage is physically 
impossible, and one must always take care not to read an evolutionary cause 
into a phenomenon that is inevitable. At the same time, it is worth posing 
the question, whether there is an optimal level of asymmetry in the segre- 
gation of damage which is non-zero. If this is true in the models, it suggests 
possibilities for future experiments, to determine whether the asymmetry is 
being actively driven by the cell, or whether possibilities for reducing the 
asymmetry are being neglected. 

In principle, this is a perfectly generic evolutionary phenomenon. It de- 
pends upon the "old-pole" experiments only to the extent that they reveal at 
least one mechanism producing such an asymmetry. Our model is intended 
to represent a large population of fissioning organisms, each of which has 
an individual "damage level" , which determines its rate of growth and divi- 
sion, and its likelihood of dying. When a cell divides, its daughters divide 
the parent's damage unequally, with the inequality controlled by a tunable 
parameter. The evolution of this population is described formally by a 
mathematical object called a "superprocess" . (A somewhat different math- 
ematical approach to the growth of structured populations was presented 
by [18] . One attempt to model the costs - but not the possible benefits - of 
asymmetries like those found in the Stewart et al. experiments is |43j.) 

An illustration of the underlying branching model is in Figure [U In this 
example, each individual cell accumulates damage at a constant rate during 
its lifetime, so that the lifelines of cells are parallel sloping lines. The only 
way to reduce the damage load is to pass the damage on disproportionately 
to one of the daughters at a division. Note the following: 

• Any fixed line of descent eventually dies. 

• The population as a whole continues to grow. 

• Death is more likely, and fissioning comes more slowly, for cells that 
are further right — that is, with more damage. 



DAMAGE SEGREGATION 
• The population eventually clusters at low levels of damage. 




► 

Increasing Damage 



FIGURE 1. Illustration of the branching model. Sloping lines 
are lifelines of cells, moving downward in time and rightward 
in damage space. Filled circles represent fission events. Open 
circles represent deaths. Arrows represent cells that survive 
to the end of the experiment. 
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At this point, we should acknowledge [75|, which appeared while the 
present paper was under review. This paper specifically models the old-pole 
effect, putting forward an account of the possible advantages of asymmetric 
damage division, which complements our more abstract model in several 
important respects. Whereas the present paper offers analytical solutions 
for models that incorporate a broad range of mechanisms for damage seg- 
regation and the effect of damage on death and splitting rates, the specific 
parametric model of [75] is amenable to simulation and computational ex- 
ploration of its particular parameter space. That model is free of the scaling 
assumptions, on which the analysis given here depends. (Roughly speak- 
ing, the scaling assumptions are analogous to those that are required for 
the validity of the diffusion approximations introduced by researchers such 
as Kimura into population genetics.) Finally, [75] differs from the present 
paper in its explicit representation of cell growth. The important lesson is 
that the same fundamental behavior is seen in two very different models, 
analyzed with very different methods. In each case, increasing the inequal- 
ity in transmission of damage to the daughter cells is found to increase the 
population growth rate in some circumstances. 

3. Description of the model 

Our model for the growth of a population of fissioning organisms such as 
E. coli is an infinite population measure- valued diffusion limit of a sequence 
of finite population branching models. Mathematical derivations of impor- 
tant properties of this model are left for section [7J and Appendix A. In the 
present section, we offer a more intuitive description of the model. 

3.1. Diffusion. The state of the process at any time t is a collection of 
cells, each of which has some level of damage. We represent the level of 
damage as a positive real number. Each cell develops independently of its 
fellows, performing four different behaviors: Damage creation, damage re- 
pair, death, and fissioning. Damage creation and repair sum to a net damage 
process, which we model as a diffusion. Diffusion processes, which are most 
commonly applied in population genetics as models for fluctuating propor- 
tions of alleles in a population (cf. [30]), may be thought of as continuous 
analogues of random walks. A real- valued diffusion is a general model of a 
random process that changes continuously in time and satisfies the Markov 
property: future behavior depends only on the current state, not on the more 
distant past. The diffusion model of damage accumulation in a cell is deter- 
mined by two parameters: the diffusion rate <r mo t(x), giving the intensity of 
random fluctuation of the damage level as a function of the current damage 
level x; and the drift b(x), giving the trend in damage, whether increasing 
(positive) or decreasing (negative), as a function of current damage. 

3.2. Branching diffusions and measure-valued diffusions. While cells 
move independently through damage space, they are also splitting in two 
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and dying, at rates that depend on their current damage state. The obvious 
way to represent the state of this "branching diffusion" process at any given 
time is as a list of a changing number of cells, each labeled with its damage 
state. It turns out to be mathematically more convenient to invert this de- 
scription and present all the damage states, each with a number (possibly 
zero) representing the number of cells in that state. Formally, such a de- 
scription is called an integer-valued measure on the space of damage states 
(that is, on the positive real numbers). A stochastic process whose state 
at any given time is one of these measures is a measure-valued (stochastic) 
process. 

We take this mathematical simplification one step further, at the expense 
of increasing the amount of hidden mathematical machinery. The discrete 
numbers of individuals in the population prevent us from applying the pow- 
erful mathematical tools of analysis. This is exactly analogous to the diffi- 
culties that arise in analyzing the long-term behavior of discrete branching 
models, or inheritance models such as the Wright-Fisher and Moran models 
[30j . The famous solution to that problem was W. Feller's [31] diffusion 
approximation. By an essentially analogous method (see, for example, [25] 
or [60]). rescaling the time and giving each cell "weight" 1/n, while letting 
the population size and the birth and death rates grow like n, and now 
sending n to infinity, we obtain in the limit a stochastic process that has 
as its state space continuous distributions of "population" over the damage 
states. Such general distributions are called measures and may be identified 
in our setting with density functions on the positive real numbers. We stress 
that the limiting stochastic process is not deterministic and can be thought 
of as an inhomogeneous cloud of mass spread over the positive real numbers 
that evolves randomly and continuously in time. Such a stochastic process is 
called a measure-valued diffusion or superprocess. An accessible introduction 
to measure- valued process produced by this sort of scaling limit, including 
the original superprocess, called "super-Brownian motion," is |66j . 

As described in [64], branching diffusions [56} [13| [29] and the related 
stepping-stone models |46[ |47| were well established in the 1960s as models 
for the geographic dispersion, mutation, and selection of rare alleles. Im- 
portant early applications of spatially structured branching models include 
[481 [65] , which helped to elucidate the spatial population distributions that 
are likely to arise from from populations under local dispersion and local 
control. The utility of measure-valued stochastic processes in population 
biology was established by the celebrated limit version of allele frequency 
evolution in a sampling-replacement model due to W. Fleming and M. Viot 
[33] , which followed quickly upon the mathematical innovations of D. Daw- 
son [13]. In particular, studies such as [271 [38] show how the measured- valued 
diffusion approach helps in taming the complexities of infinite sites models. 
Superprocesses have found broad application in the study of persistence 
and extinction properties of populations evolving in physical space (such as 
[78] [26] ) or in metaphorical spaces of genotypes (such as [191 E3 CEZI ) • A 
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thorough introduction to the subject, including extensive notes on its roots 
and applications in branching-diffusion and sampling-replacement models, 
is US]. 

3.3. The damage-state superprocess. Formally, we posit a sequence of 
branching models, converging to our final superprocess model. If the scaling 
assumptions are reasonable, important characteristics of the branching mod- 
els — in particular, the connection between the process parameters and the 
long-term growth rate — should be well approximated by the corresponding 
characteristics of the superprocess. 

The ra th model in the sequence begins with a possibly random number 
N n of individuals. Individuals are located at points of the positive real 
line ffi + . The location of an individual is a measure of the individual's 
degree of cellular damage. We code the ensemble of locations at time t > 
as a measure X™ on IR + by placing mass 1/n at each location, so that 
the initial disposition of masses Xq is a random measure with total mass 
N n /n. During its life span the damage level fluctuates, as the organism 
accumulates and repairs damage. An individual's level of damage evolves as 
an independent diffusion on ]R + with continuous spatially varying drift b(x), 
continuous diffusion rate cr mot (x), and continuous killing rate (depending 
on n) k n (x). While our mathematical framework allows general behavior 
(mixtures of reflection and killing) at the boundary point 0, it seems sensible 
for applications to assume complete reflection. Other boundary behavior 
would imply a singular mortality mechanism at damage level 0. 

An individual can split into two descendants, which happens at spatially 
varying rate a n (x). As in the standard Dawson- Watanabe superprocess 
construction [28j [16j [25] , we assume that birth and death rates grow with n 
but are nearly balanced, so that there is an asymptotically finite net birth 
rate 

(1) (3{x) := lim (a n (x) - k n (x)) < oo 

n— >oo 

while the total rate at at which either death or splitting occurs satisfies 

(2) p(x) := lim n~ l (a n (x) + k n (x)) < oo. 

n— »oo 

In other words, an individual at location x attempts to split at a rate which is 
asymptotically np(x), and such an attempt is successful (resp. unsuccessful) 
with a probability that is asymptotically \(l+(3{x) /n) (resp. ^(1— (3{x)/n)). 
We assume that (3* := sup x /3(x) < co and fioo := lim^oo f3{x) exists 
(though it may be — oo). Both of these conditions hold under the biolog- 
ically reasonable assumption that (3 is non-increasing (that is, asymptotic 
net birth rate is non- increasing in the amount of damage). 

The unbalanced transmission of damaged material leads to one important 
difference from the standard Dawson- Watanabe superprocess construction: 
When a split occurs the two descendants are not born at the same point 
as the parent. A parent at location x has descendants at locations x — Y 
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and x + Y, where Y is randomly chosen from a distribution 7r" on K + . We 
suppose that 

/oo 
y 2 dir^iy) < oo 
-oo 

for some continuous function a 2 eg . The quantity o- 2 cg (x) is a measure of the 
amount of damage segregation that takes place when a parent with level of 
damage x splits. 

We show in the Appendix that the process X n converges to a limit process 
X, as n — > oo, as long as the initial damage distributions Xq converge to 
some finite measure Xq. The limit is a Dawson- Watanabe superprocess 
whose underlying spatial motion is a diffusion with drift b (the same as the 
original motion) and diffusion rate a, where 

(4) a 2 {x) := a^ ot {x) + cr 2 cg (x)p(x) 

Thus the asymptotic effect of segregation is equivalent to an increase in the 
diffusion rate, which may be thought of as the extent of random fluctuation 
in the the underlying motion of an individual cell through damage space. 
The same effect may be achieved by increasing simultaneously the fission 
and death rates, leaving fixed the net rate of birth. The quantity o- 2 cg (x) 
measures a re-scaled per generation degree of damage segregation, and multi- 
plication by the re-scaled splitting rate p(x) produces a quantity a 2 cg (x)p(x) 
that measures a re-scaled per unit time rate of damage segregation. 



4. Major results 

Write Z± for the asymptotic re-scaled population size at time t. That is, 
Z t is the total mass of the measure X t . We show that there is an asymptotic 
population growth rate X in the sense that the asymptotic behavior of the 
expectation of Zt satisfies 

(5) lim t -1 bgE[Z t l = A. 

t~+oo 

Thus, E[Z t ] grows to first order like e— . 

When the asymptotic growth rate A is positive, this result may be sig- 
nificantly strengthened. Theorem 17.41 tells us then (under mild technical 
conditions) that: 

• To first order, the total population size at time t, Zt, grows like e-'. 
More precisely, Zt/E[Zt] converges to a finite, non-zero, random 
limit as t — > oo. 

• The random distribution describing the relative proportions of dam- 
age states of cells at time t converges to a non-random distribution 
with density ccj) as t — > oo, where <f> is the unique solution to the 
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ordinary differential equation 

\ [a 2 (x)<f>(x)] " - [b(x)cf>(x)}' + (3(x)<f>(x) = M 

with ^V(O) = 26(0)0(0) 

and c is the normalization constant required to make c(p integrate 1. 
That is, for any < a < b < oo the proportion of the population with 
damage states between a and b converges to f <j)(x) dxj J °° <p(x) dx. 

Suppose now that we identify (as is commonly done) fitness of a cell type 
with the asymptotic growth rate A. It is not obvious, under arbitrary fixed 
conditions of average damage accumulation and repair (6), splitting rate (p), 
and net reproduction (/3), what level of randomness in damage accumula- 
tion and repair (cr mo t) and damage segregation (cr seg ) has the highest fitness. 
Under the biologically reasonable assumption that the average damage ac- 
cumulation rate is everywhere positive, though, it is generally true that 
moderate degrees of randomness in both mechanisms lead to higher growth 
rates than either very small or very large degrees of randomness. Theorem 
17.51 (in section I7.3P tells us that, under fairly general conditions, if the dif- 
fusion rate a is reduced uniformly to 0, or increased uniformly to oo, this 
pushes the asymptotic growth rate down to = lim z _ i . tX3 (3(z). 

Since the growth rate at any damage level is higher than flao, this result 
tells us that the growth rate is minimized by sending a uniformly either to 
or to oo. Too much or too little randomness in damage accumulation yields 
lower growth rates; optimal growth must be found at intermediate levels of 
diffusion. 

Imagine now a situation in which the upward drift in damage has been 
brought as low as possible by natural selection, optimizing repair processes 
and tamping down damage generation. We see that there may still be some 
latitude for natural selection to increase the lineage's growth rate, by tin- 
kering with this randomness. Recall that a is composed of two parts: the 
randomness of damage accumulation during the lifecourse, and the inequal- 
ity of damage transmission to the daughters. Without altering damage gen- 
eration or repair, but only the equality of transmission to the offspring, the 
cell type's long-term growth rate may be increased. In particular, if dam- 
age accumulation rates are positive and deterministic — a fixed rate b(x) of 
damage accumulation per unit of time, when a cell has damage level x - 
the growth rate will always be improved by unequal damage transmission. 

5. An example 

Suppose that damage accumulates deterministically during the life of an 
organism, at a rate proportional to the current level of damage. In the 
absence of damage segregation, then, the damage at time t grows like e bt 
(this corresponds to b(x) = bx and c^ot = 0). We suppose that the combined 
parameter a 2 p{x) describing the time rate of damage segregation for an 
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organism that fissions at damage level x is r 2 x 2 , where r is a tunable non 



negative parameter (hence, a(x) = ya 2 cg p(x) = tx). The generator for the 
spatial motion in the limit superprocess is thus 

f. x 2f»( x ) +bx f( x ). 

This is the generator of the geometric Brownian motion exp(rBt + (b — 

2 

We see immediately that adding enough noise can push the spatial 
motion toward regions with less damage. To make this fit formally into our 
theorems, we assume complete reflection at 1. (We have previously taken 
the range of damage states to be the positive real numbers (0, oo), whereas 
here we have taken the interval (l,oo). The choice of interval is irrelevant 
as far as the mathematics goes, and we could convert this choice to our 
previous choice by simply re-defining b{x) and cr(x)). 

Suppose that the bias of reproduction over dying in the superprocess is 
of the form (3{x) = (3$ — (3\x. The population growth rate will be the same 
as for the process Yt = Bt + — ^)t, with reflection at 0, and killing at rate 
/?o — Pie rYt ■ The eigenvalue equation is 

+ bxf'(y) + (A> - 0ie Ty )f(y) = A/(y), 

where b := ^ — 77, with boundary condition /'(0) = 0. 

Using standard transformations as in section 6 of [67], or symbolic math- 
ematics software such as Maple, we may compute a basis for the space of 
solutions to be the modified Bessel functions 



Iu (v/8/3i/r 2 e ry/2 ) e~' by and K u (y8/3iA 2 e ry/2 ) e 



by 



where 



(6) 



2JP - 2(A) - A) 



By definition, A is the infimum of those A such that the solution / is 
non-negative. Since /'(0) = 0, it cannot be the case that /(0) = 0. Hence, 
the solution for A has no zero, but the solution for every A < A does have 
a zero in M + . By continuity (as argued for the same problem in [67]), it 
follows that the solution for A is a multiple of K v . The boundary condition 
translates to 

( 7) (vm\ = _L 

{ ' K v \ t ) 

Find uq real such that the last zero of Ki UQ is at i/8/3i/r, and represent 
the corresponding A by 

u 8 2 MU 
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Damage segregation parameter (t) 

Figure 2. Asymptotic growth rates A computed from (jlip 
for populations whose deterministic rate of damage accumu- 
lation is equal to the current damage level, and whose relative 
excess of births over deaths is 1.9 — x/2 at damage level x. 
Damage segregation is added at rate tx for births to parents 
at damage level x. 



Then 

As an example, take (3\ = |, b = 1 and /% = 1.9. Using the mathematical 
programming language Maple we compute the asymptotic growth rates given 
in Figure [21 The maximum growth rate of occurs at r = 1.75. Figure [3] shows 
the asymptotic distribution of damage for different values of r. Note that 
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Figure 3. Asymptotic distribution of damage for the choice 
of parameters j3\ = \,b = 1 and (3q = 1.9, and various choices 
of r. 



as r increases, the asymptotic distribution becomes increasingly compressed 
into the lower range of damages. 

Remember that a = J cr^ot + a scgPi where cr^ot quantifies the random- 
ness of fluctuations in the damage during the life of a cell; a^ cg quantifies 
the extent to which the damage is split unequally between the daughters 
and p is the splitting rate. More unequal distribution of damage among the 
daughter cells means that there will be effectively more randomness in the 
rate of damage accumulation within a line. 

We point out one implication for the evolution of generation time. Sup- 
pose that there is no mechanism available for altering the extent of damage 
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segregation in each fission event. (Perhaps damage segregation is not ac- 
tively determined, but is only a consequence of the inevitable one-way in- 
heritance of an "old pole", or simple random fluctuation.) Suppose, though, 
that the organism can accelerate the generations. This has the effect of in- 
creasing a, which may be enough to shift the asymptotic growth rate into 
the positive. 

To be specific, suppose there is a population of 1000 bacteria, with de- 
terministic damage accumulation at a rate x units per day when the cell 
has x units of damage. A cell with x units of damage has a hazard rate for 
reproduction in the next instant of 25.95 — x/4, and a hazard rate for dying 
of 24.05 + x/4. (Measuring time in days, this corresponds to a mean time 
between fission events of about 1 hour.) In the notation of our model, this 
means that p{x) = 50/1000 = 0.05 and f3(x) = 1.9 — x/2. The minimum 
level of damage in a cell is 1 unit. At a fission event, if the mother has 
x units of damage, one daughter gets 90% of x and the other 110% of x. 
(Damage units must be thought of as proportional to the total volume of the 
cell. Thus, when a cell with x damage units divides, there is 2x damage to 
be divided up.) From equation ([3]), we compute that <7 seg (x) = 3.1x. Since 
Cmot was assumed to be 0, this means that a(x) = a scs (x)\/~X)5 = 0.7x. We 
see from Figure [2] that the asymptotic growth rate of the population will be 
—0.09, meaning that the population may be expected to decline over time. 

Suppose now there were a mutant strain which had boosted its rate of 
reproduction to 40 — x/4, but was penalized with a rate of dying 38.2 + x/4. 
That is, fission rates were increased by 14.05 (about 35% at the minimum 
damage state), but the mortality rate was increased by a bit more, by 14.15. 
Everything else stays the same. We now have p = 78.2/1000 = .0782, 
meaning that cr(x) = 0.87x, and despite the mortality penalty we see that 
the population now has a positive asymptotic growth rate of 0.09. Instead 
of declining, this strain will tend to double in numbers about every 10 days. 

6. Conclusions 

We have introduced a mathematical model of population growth, for a 
population of haploid organisms that accumulate damage. This is a limit 
of branching processes, in which individuals accumulate damage, rising and 
falling at random, with an upward tendency. Under fairly general assump- 
tions we have shown 

• The population growth rate converges to a fixed rate, determined by 
the solutions to an ordinary differential equation (ODE); 

• If the asymptotic growth rate is positive, the relative proportions of 
damage levels within the population converge to a fixed distribution, 
also given by the solution to an ODE; 

• The effect of increasing damage segregation in the model (that is, 
causing one daughter cell to have higher damage than the parent 
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after the split and one lower) is equivalent to increasing the random- 
ness (diffusion) of the damage-accumulation process; 

• Accelerating the turnover of generations — increasing birth and 
death rates equally, leaving the net birth rate unchanged — effec- 
tively increases damage segregation; 

• In general, the optimum level of combined damage diffusion and 
damage segregation (as measured by the growth rate) is not 0, but 
is some finite non-zero level. 

This is one possible mathematical analysis of exogenous repair. A strain 
of cells with deterministic damage accumulation will eventually run itself to 
extinction because of the increasing damage load. Perhaps surprisingly, the 
simple expedient of dividing the damage unequally between the offspring 
may be enough to rescue the line. In general, if the damage diffusion rate 
were sub-optimal, a mutant line could obtain a selective advantage by in- 
creasing the level of damage segregation. In fact, damage segregation could 
be a worthwhile investment even if it came at the cost of a slight overall 
increase in damage accumulation. Perhaps even more surprising, simply 
accelerating cell division and cell death equally, leaving the net birth rate 
unchanged, may be enough to shift a negative population growth rate into 
the positive range. Of course, there can also be too much damage segre- 
gation, and there is potential for improving growth rates by equalizing the 
inheritance between the daughter cells, if the variation in damage accumu- 
lation is excessive. 

An abstract model such as presented in this paper is primarily a guide to 
thinking about experimental results, rather than a template for analyzing 
the data from any given experiment. This does not mean that this model 
cannot be experimentally tested. Some general predictions of the model are: 

(1) Mechanisms for actively increasing damage segregation in reproduc- 
tion should be common. 

(2) While individual lines tend toward increasing damage, the overall 
distribution of damage in the population stabilizes. 

(3) Damage segregation compensates for "too little" randomness in 
damage accumulation. Thus more unequal division of damage might 
be expected in organisms whose rate of damage accumulation is more 
deterministic. 

(4) Manipulating the extent of damage segregation should affect the 
population growth rate. Under some circumstances, increasing seg- 
regation should increase the growth rate. 

(5) One way of manipulating damage segregation is to manipulate the 
generation time. Increasing the fission rate and the death rate 
equally has the effect of increasing the segregation effect. 

The first "prediction" is already fairly established. While asymmetric 
reproduction and exogenous repair were implicit in the century-old line of 
experiments recounted in [8], the direct measure of asymmetry in unicellular 
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reproduction has only recently become possible. Some of the organisms 
studied are 

(1) Caulobacter cresentus: This bacterium grows from one end and then 
divides. The old end accumulates damage and senesces, while the 
new end experiences "rejuvenating reproduction" [3]. 

(2) Escherichia coli: While this bacterium seems to divide symmetri- 
cally, rejuvenation seems to happen preferentially in the middle, 
while the poles accumulate damage. [70] found declining vigor with 
increasing age of the inherited cellular poles. 

(3) Saccharomyces cerevisiae: The aging of the mother cell in these 
budding yeast was established as long ago as [57], and the func- 
tional asymmetry between mother and daughter was established by 
[41j . [51] showed that the continuing lineage growth depends upon 
maintaining the asymmetry in damage state between mother (more 
damaged) and daughter (less damaged). An atp2 mutant strain that 
failed to segregate active mitochondria preferentially to the daugh- 
ter cell succumbed to clonal senescence. [I] found that carbonylated 
proteins are selectively retained by the mother cell, and that this 
requires an active SIR2 gene. (Higher activity of SIR2 also increases 
the longevity of the mother cell itself, which seems surprising on the 
face of it; but this is likely a consequence of epistatic effects, cf. [44].) 

(4) Schizo saccharomyces pombe: These fissioning yeast seem superfi- 
cially to divide symmetrically. While there is no gross morphological 
distinction on the basis of which the two products of reproduction 
could be allocated to categories "parent" and "offspring," the fis- 
sion, there is heritable asymmetry between the two fission products 
in size [7] and number of fission scars [10] . The inheritance of dam- 
aged proteins was studied in [55] , but asymmetry of inheritance was 
not measured directly. 

Yeast seem to be a promising class of organisms for studying the inheri- 
tance of damage. As remarked by Lai et al. [51] . 

The aging process of yeast cells best illustrates the cellu- 
lar and generational asymmetry [....] The immortality of 
the yeast clone] depends upon the establishment of an age 
difference between a daughter cell and the mother cell that 
produces it. Although the daughter cell receives cellular com- 
ponents from the mother cell, it does not inherit the mother 
cells age and is always younger than the mother cell. This 
age asymmetry implies that any substance or process respon- 
sible for the aging of mother cells must be carefully isolated 
from the daughter cells. 

The biodemography of both S. cerevisiae and S. pombe is relatively well un- 
derstood, and methods for tracking fission scars, carbonylated proteins and 
inactive mitochondria are well developed. We have two species with quite 
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different types of damage segregation, and comparisons are likely to prove 
illuminating. The identification of other clonal senescence mutants that act 
on damage segregation, whether variant mutations to SIR2 or another gene, 
would expand the toolkit for tuning damage segregation in S. pombe. 

The recent study of carbonylation [55] seems particularly amenable to 
analysis in terms of exogenous repair and our model. This study examined 
the distribution of damaged proteins in the whole population, finding that 
the dividing cells were concentrated down near 0, while the non-dividing 
cells had a significantly higher level of damage. The main parameters in our 
model — the rate of damage accumulation in stationary phase, the vital rates 
for cells with a given damage level, the type of damage inheritance — could 
all be measured in principle. Indeed, the symmetry of damage inheritance 
between the daughters was examined, but only to point out that there is 
some sharing of damage, that the carbonylated proteins are not all retained 
by one of the two daughters. Pushing this further, it would be possible to 
measure the inequality of damage inheritance, both of carbonylated proteins 
and inactive mitochondria. Once these parameters had been estimated, it 
would be possible to compare the damage distribution in the paper with our 
theoretically predicted asymptotic distribution. 

Comparing mutant strains with different damage segregation properties 
would further refine the comparison. It would also be possible to perform 
competition experiments: Broadly speaking, we would predict that as the 
overall damage load increases — perhaps by exposure to paraquat or dis- 
abling of antioxidants — this would shift the optimal level of damage segre- 
gation. In some simple choices of parameters for the model this shift would 
be toward more extreme damage segregation, the same effect that would be 
predicted from decreasing the variability of damage accumulation. Thus, 
we might expect that strains exposed to such toxic environments over some 
time would evolve to show higher levels of damage segregation. If damage 
segregation is not controlled by active mechanisms, this means that increas- 
ing damage accumulation should select for a faster turnover of generations: 
simultaneously higher rates of cell death and fissioning. 

Perhaps no simple resolution should be expected. We are examining the 
interaction of individual-level aging downward, with aging of parts of indi- 
viduals, and simultaneously upward, toward population-level aging. In the 
present model one might say that there is a partial individualization of senes- 
cence. That is, while there is no identifiable individual who ages and dies, 
leaving behind youthful descendants, damage segregation has transformed 
the senescence problem plaguing the entire population into a problem of 
some individual cells. 

There has been some controversy (see [77], and the response [71]) about 
whether results such as those of |70j reflect genuine "aging." Perhaps more 
profitable is to see, in this growing family of experimental results, clues to the 
broader context of aging: Paradigmatic aging in our species and similar ones 
is one of a class of mechanisms which function to redistribute damage within 
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a population. Most metazoans, perhaps inevitably, have converged on the 
extreme strategy of keeping essentially all the damage (even generating much 
more damage in the process), and produce pristine offspring. In unicellular 
organisms, there seem to be a broader range of mechanisms, and degrees of 
damage differentiation in the end-products of reproduction. 

Is this relativization of senescence of any relevance for metazoans, then? 
It is, for at least two separate reasons. The first is that the propagation of 
the germ line is in many senses comparable to the propagation of a line of 
unicellular organisms. As eloquently commented by Lai et al. [51], "Some 
of the key questions in aging concern the differences between germline and 
soma. Any mechanism invoked to explain the aging of the soma must also 
be able to accommodate the immortality of the germline. At the level of 
the organism, the issue is the renewal that occurs at each generation, pro- 
viding the progeny a fresh start with the potential for a full life span." This 
principle may also be significant for the regulation of replicative senescence 
in somatic stem cells or other somatic cells. The recent discovery of non- 
random segregation in mouse mitotic chromosome replication [5j [6] may be 
a hint in this direction. Even more suggestive is the apparent functioning 
of this mechanism in maintaining the integrity of cell lines during the lives 
of individual higher eukaryotes (small intestine epithelial crypt cells in H. 
sapiens and neural precursor cells in D. melanog aster) [63| . 

The other reason is perhaps more profound, and certainly more specula- 
tive. Over the past several decades, mathematical theory has played a sig- 
nificant part in the evolving discussion of risk-management and bet-hedging 
in natural selection. J. Gillespie's seminal treatment [36] notes that 

the variance in the numbers of offspring of a genotype has 
two components, the within-generation component resulting 
from different individuals of the same genotype having dif- 
ferent numbers of offspring, and a between-generation com- 
ponent due to the effects of a changing environment. 

The general conclusion has been that lower variance in essential demographic 
traits is selectively favored [341 [37] , although this picture is complicated by 
varying environments, if the variance is associated with variable response to 
the environment [22J. Damage accumulation falls into a third class of vari- 
ance: variable epigenetic inheritance of a resource shared between offspring. 

Our model departs from the models of within-generation variation in 
which the varying phenotypes of siblings are independent, but is allied 
with models of resource sharing within families. In the case of damage- 
accumulation the "resource" is the undamaged cell components, but this 
model might in principle be applied to the transmission of resources or sta- 
tus from parents to offspring. In contrast to the general selective advantage 
for reduced within-generation variability in the uncorrelated setting, the re- 
sults of our analysis suggest that increased variability may be selectively 
favored when the phenotype is a shared resource. If the resource is simply 
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size, then our model might be compared to the asymmetric-division cell cycle 
model developed by [73]. That work lacked an explicit stochastic popula- 
tion model, and did not address the question of overall population growth 
rate, but some of the observations from that paper may be relevant, most 
significantly the simple empirical fact of size asymmetry in yeast division. 

This beneficial variance may have some relevance for the growing recog- 
nition that the evolution of longevity depends fundamentally on the nature 
of intergenerational transfers. All the work on this problem of which we 
are aware, particularly [52l EU [TT] , implicitly assumes that resources are 
divided equally according to need (though [12] points out that the goal of 
equality may still produce systematic biases based on birth-order). The key 
paper of Trivers [72] explicitly opposes a parental goal of equality against 
the offspring's interest in monopolizing resources. Models of offspring pro- 
visioning in this tradition, such as [58^ I62|. seem to assume a "diminishing 
returns" to investment, hence higher fitness for reduced spread in offspring 
quality, except when there is a threshold for survival that only one is likely to 
cross. Our model suggests that the cumulative effect of exponential growth 
can easily make unequal division of resources a winning strategy. This in- 
equality may be purely random, and need not depend on sibling competition 
or some other strategy for identifying and rewarding the inherently superior 
offspring. 

We note that intergenerational transfer theory must address the com- 
plication of parent-offspring and sibling competition, because of the genetic 
differences between relatives. Our model avoids those complications, assum- 
ing a complete comity of cells. The parent gives everything to its offspring, 
and the daughters are as interested in the other's welfare as in their own. 
Nonetheless — and this is what may be surprising — equal endowment of 
offspring may not be optimal. A daughter's fitness may be best served by 
trading the sure 50% inheritance for a lottery ticket that may yield more, 
although it may yield less. 

7. Mathematical methods 
Define (Pt)t>o to be the semigroup generated by 

£0 := + b<f>' + {3<p 

with boundary condition 

(8) Po^^'(O) = (l-Po)0(O) 

for some fixed < po < 1. (This semigroup will not, in general, be sub- 
Mar kovian because (3 can take positive values). As we observe in Section 
17.21 K[(<fi, X(t))] = L + Pt4>(x) dv{x) when Xq is the deterministic measure 
z/, and so the behavior of {X t )t>o is governed by that of (Pt)t>o, at least at 
the level of expectations. The choice po = 1 is most relevant for the descrip- 
tion of fissioning organisms, because it corresponds to complete reflection 
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at the state representing no damage. Other choices of po would imply the 
rather anomalous presence of an additional singular killing mechanism at 
0. However, since the mathematical development is unchanged by the as- 
sumption of a general boundary condition, we do not specialize to complete 
reflection. 

It turns out that the entire long-term behavior of (Xt)t>o is rather simply 
connected to that of (Pt)t>o- Let (P^)t>o be the adjoint semigroup generated 
by 

with boundary condition 

(9) (0'(O) - 26(0)0(0)) = (1 - po)0(O). 

For any A, define 4>\ to be the unique solution to the initial value problem 
(10) 

L*<j> x = with <j> x (0) = Po and £^L(^(0) - 26(0)^(0)) = 1 - p . 

Put 

(11) A := infjA : (f>\(x) is non-negative for all x\. 

We recall in Proposition 17. 21 that, under suitable hypotheses, for a compactly 
supported finite measure v 

(P t+S l,v) _„ As 



lim — — — = e 

f^oo (P t l,u) 



so that a fortiori 
and, moreover, 



lim t _1 log(P t l, v) = A, 

t— »oo 

lim v ) = /o°° <-P(z)<P\(z)dz 

£2o{P t l,v)~ !™<t>\{z)dz 



10 

for a bounded test function (p. We subsequently show in Theorem 17.41 again 
under suitable hypotheses, that if A > 0, then 

(<p,X t ) <p(z)<l>x(z)dz 

t^oo(P t l,u) f °°Mz)dz 

in probability, where W is a non-trivial random variable (that doesn't de- 
pend on ip). In particular, the asymptotic growth rate of the total mass 
(1, Xt) is A and 

lim (gijfj) = /o°° <P( z )<l>\( z ) dz 
t™{l,X t ) f °°Mz)dz 

We show in Section 17.31 under quite general conditions (including the as- 
sumption that (5 is non- increasing so that (3^ = inf x f3(x)) that the presence 
of some randomness in either the damage and accumulation diffusion or the 
damage segregation mechanism is beneficial for the long-term growth of the 
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population but too much randomness is counterproductive. That is, if b, (3 
and p are held fixed, then A converges to Poo as the diffusion rate function 

a = J ci^ lot + (Xg Cg p goes to either or oo, whereas A > [3^ for finite a. 

7.1. Results from [67]. Essential to linking the long-term growth be- 
haviour of the superprocess to the theory of quasistationary distributions 
is the observation that L — (3* (recall that j3* is a constant, the upper bound 
on (3) is the generator of a killed diffusion semigroup Qt4> = t Pt4>- Thus, 
when we have a result stating that Qt converges to a stationary distribution 
with density a multiple of cp\, meaning that for any bounded test function 
(p : ]R + , and any compactly supported finite measure v, 

t->cc (Q t l,u) f£°4>x(z)dz 

then the same holds for Ft- 

For t > and v a compactly supported finite measure on K + , set 

(12) ^):- {PM 



m(l) • 

When v = 5 X , write u>t(x). Note that 



(13) u t (u) = / uj t (x) du(x 







for any compactly supported finite measure v. Define 

(14) tu*(x,v) = sup — rT . 

t > uJ t (y) 

Assume the following condition: For all compactly supported finite mea- 
sures v on M + , and some e > 

(GB2) sup sup e" (2A " e)s / uj*{y,v) 2 P s (x,dy) < oo. 

s>0 xSK+ J 

While (GB2) is an abstract condition, it is implied by a fairly concrete 
condition on the diffusion parameters. 

Lemma 7.1. Condition (GB2) always holds if the net birth rate (3 is mono- 
tone non-increasing, and there is complete reflection at 0. 

Proof. Changing (3 by a constant leaves u>t, hence u* as well, unchanged. 
Thus, we may assume that (3 < 0, so that Pt is the semigroup of a killed 
diffusion. 

Given x < x', we may couple the diffusion Xt started at x with a diffusion 
X' t , with identical dynamics but started at x' , so that Xt < X[ for all 
t. By the monotonicity of f3, this implies that Ptl(x) > P t l(x'), since 
PtX{x) = F x {tq > t}. Consequently, uj t is monotonically non- increasing, 
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and so is w*. In particular, for all y, u*(y,v) < u*(Q,v), which is finite, by 
the results of section 2.5 of [67] . Then 



sup sup e ( 2 - ^ s / a>*(y, v) 2 P s {x, dy) < w*(0, v) 2 < oo. 

□ 

We will also need to assume the technical condition 
(LP) £j* is in the limit-point case at oo. 

We define 

~ boF- l {x) 1 , x 

b{z) := aoF-\ X ) -r° F 

and (3(x) := j3 o F~ l (x), where F(y) := du/o~{u) is the Liouville trans- 
form, the spatial transformation that converts the diffusion coefficient into 
a constant. Then Lemma 2.1 of [67] says that (LP) holds whenever 

(LP') lim inf x~ 2 (b{x) 2 + b'{x) + 2/?(x)) > -oo. 

Examples which do not satisfy this condition must have pathological fluctu- 
ations in b. 

Define ip\ to be the unique solution to the initial value problem 

(15) Lipx = \?px with Va(0) = Po and ^-^(Q) = 1-po- 

We summarize relevant results from Lemma 2.2, Theorem 3.4 and Lemma 
5.2 of E3 



Proposition 7.2. Assume the conditions (GB) and (LP). 

(i) The eigenvalue A is finite. 

(ii) The semigroup (Pt)t>o has an asymptotic growth rate r\. That is, for 
any compactly supported v, and any positive s, 

(16) lim = e <». 

V ^ (P t l,u) 

(iii) The following implications hold: 

A > Ax> V = A > (3oo 

A = /3oo V = A = Poo 

A < /?oo => f? = A or 7] = /3oo 

(iv) // r/ = A ^ /?oo, i/ien, /or any compactly supported finite measure v, 

(17) ]imu*(„)= / ^TT^(^), 

i^oo 7 ^(1) 

and, for any bounded test function (p, 

(18) Um ( p tV, v) _ / °° <p(z)<j>x{z)dz 



t^(P t l,v) !™ct>x{z)dz 
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(v) If (3 is non-constant and non-increasing and po = 1 (that is, there 
is complete reflection at ), then (|18h holds if and only if A > /3oo • 
In particular, f|18j) holds if po = 1 and /? is non-increasing and un- 
bounded below (so that (3^ = —oo). 

Proof. The only part that is not copied directly from |67| is (v) . If A > f3oo , 
then n = A > (3^ by (iii) and (fT8j) holds by (iv). Suppose that (fTKl) holds but 
A < (3<x>- Since there is no killing at (po = 1), by ([8]), the constant function 
is in the domain of XL. Thus, we may apply the Kolmogorov backward 
equation to see that by ([TBI) . 



- 'L 
^ ds 

- <L 
ds 

- *L 
ds 



\t~*co (P t l, u) 



lim 

s=0 \t— »oo 



J °° P a l(z)(f>x(z)dz 

V 0A,(z)dz 



by CI 



_ / °°x:i(z)(/>A(z)dz 

Jo°° <Mz)<fe 
_ !™ P{z)<l>\(z)dz 
/o°° 0A(«)<fo 

Hence 77 > /3oo because /3 is non-constant and </>a(^) is non-negative for all 
z and strictly positive for z close to 0. (In fact, <p\ is strictly positive for all 
z). However, this contradicts (iii). □ 

A simple consequence of Proposition 17.21 is the following. 

Corollary 7.3. For any compactly supported finite measure v and any e > 0, 

there are positive constants < c ej „ < C €ji , < 00 such that if < t' < t" < 
00, then 



(19) 



D (r,-e)(t"-t') < (Pt"l,v) < c 



o (n+e)(t"-t>) 



Proof. By replacing P t by e~^ 3 * t Pt, we may suppose that (Pt)t>o is the semi- 
group of a killed diffusion, so that (.Ptl, v) is non- increasing in t. 
There is some integer N such that if iV < n is an integer, then 



{P n i,y) 

If N < n' < n" are integers, then, since 

n"-l 

n 



< e^ +e . 



(Pn»l,v) 
(PrA,v) 
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it follows that 

(T,-e)(n"-n') < (^n" 1 ^) < <n+e)(n" -n 1 ) 

Note that if n', n" are arbitrary non- negative integers, then 

(P n "l, v) _ {P n "l,v) (P n >s/ N l,v) {P n >i vN l,u) 
(P n 'l,v) (P n n VN l,u) (P n il,u) {P n > yN l,v)' 

On the one hand, 

(P m »l,u) {Pn1,v) , (P n "l,v) (P n , vN l,v) 
mm — — r < 



m',m"<N (P N 1,U) (P m 'l,v) (P n »wNh^) (^n' 1 ^) 

(P m »l,u) (P N 1,U 



< max 

m',m"<N (PnI 



On the other hand, 



e -N\ V -e\ e ( V -e)(n"-n') < (-Pn"Vjvl, v) < e JV|^+ £ | e (^+ £ )(„''_„') _ 

(P n t VN l,u) 

Thus the claimed inequality holds with suitable constants for t', t" restricted 
to the non-negative integers. 

For arbitrary t' < t" with \t'~\ < [t"\ we have 

c ^ e -2\ v -e\ e ( V -e)(t"-t') < H([n-[fJ) 

< (P^l,^) < (P f »l,i/) < gtfj 

- <iVji,i/> - (Pt/i.i/) - (p^M) 

and the result holds for such t',t" by suitably adjusting the constants. 
If t' < t" but \t'] > [t"\ (so that t" - t' < 1), then the observation 

(P^l,^- 1 < (P f „l,z.) < f(%jV ; 



(p Lt „ji,i/)y - (p t ,i,i/) - v^ViM 

shows that a further adjustment of the constants suffices to complete the 
proof. □ 

7.2. Scaling limit. The asymptotic behavior of Dawson- Watanabe super- 
processes has been undertaken by Englander and Turaev [23], who applied 
general spectral theory to demonstrate, under certain formal conditions, the 
convergence in distribution of the rescaled measure e~ Xct Xt, where A c is the 
generalized principal eigenvalue. Convergence in distribution was improved 
to convergence in probability by Englander and Winter 23J. We take a 
different route to proving a similar scaling limit, which has two advantages. 

(1) Our proof is more direct, and the conditions, when they hold, 
straightforward to check. 

(2) Our proof holds in some cases when the results of [24] do not apply. 
In particular, the scaling need not be exactly exponential. 
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On the other hand, our approach is more restrictive, in being applicable 
only to processes on R + . 

Our proof depends on Dynkin's formulae [2T| for the moments of a 
Dawson- Watanabe superprocess. For any bounded measurable ip, tp 1 : M + — ► 
[— 1, 1], any starting distribution u, and any < t < t', 

E v [(<p,X t )] = (P t tp,u), 
(2Q) E v [(<p,X t WM] = (P t <p,v)(P v <ff,v) 

+ (/ P s [p- P t - S <p- P t >-s<p'] ds,v). 
Jo 

Note by Proposition 17.21 that in the following result the assumptions rj = 
A > and (US) are implied by the assumptions i] — A > and A ^ (5^ . 



Theorem 7.4. Suppose that the conditions (LP) and (GB2) of Section 7.1 



hold. Suppose further that r\ = A > and (|18p holds. Then the rescaled 
random measure Xt : = converges in L 2 to a random multi- 

ple of the deterministic finite measure which has density 4>\ with respect to 
Lebesgue measure. That is, 

(21) W:=\im(l,X t ) 



exists in L 2 , and if ip is any bounded test function, 
(22) limE„ \(<p,X t )-W 



J °° ip(z)(j)x(.z) dz 



io <P\( z ) dz 

The long-term growth rate of the total mass of Xt is lim^oo t~ l log(l, Xt) 

A. 

Proof. From Proposition 17.21 we know that 

{P t <p,v) _ Jo" f(z)(f>x(z)dz 



Also, 



Km E u [((p,X t )} - 
t->oo {Pthv) Jo 4>x{z)dz 



(Pti^y^p^i^y 1 [ p s [ P -Pt- s vP t >- s ip'] {x)ds 

Jo 

(Pt- S l,v) (Pf-sl,") 



(P t l,u) (P*l,v) 



ut-sjy) ut>-a(y) Pts<p{y) Py- a <p (y) p , , x , 

x P(y) TT 7 \~B — 7TT~5 ——Ps{x,dy)ds. 

ujt-s{y)ut<-s[v) Pt-sHy) Pf-aHy) 

For fixed y the integrand converges as t, t' — > oo to 

n(„-\ ( t^l \ 2 ( fo° <P( Z )M z ) dz \ 2 -2Xs 
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By (|T9l) . for any t < t' the integrand is bounded for any positive e by 

* 1 -2s(X-e) ( u t(v)\ 2 ^ * 1 -2s(\- e ) i \2 

p e v - ' sup — — < p e v - >uj*{y, v) , 

C t ,u V r U T (V) ) C e>ll 

where p* := swp x p{x). Condition (GB2) implies that this upper bound has 
finite integral for e sufficiently small. This allows us to apply the Dominated 
Convergence Theorem to conclude that 

lim (P t l,i/)- 1 (P t ,,l,i/}- 1 //' P s [p ■ P t _ a <p ■ P t ,- S ip') ds, ; 
it->°° \Jo 



(23) 



t'>t 



J °° <p(z)<l>\(z)dz 



~2Xs 



(P s ptl>l,i>)ds. 



\$h> v ) Io° ^ z ) dz ' Jo 
That is, for any e > 0, there is T such that for all t' > t > T 

(pi, I/)" 1 <Pt/, 1, z^ 1 ^ jf P s [p • P t _ a¥ > • Pt/_ a y/] ds, - M{v, ip) 

where M{u,ip) is the right-hand side of ()23l) . 

By d20J), we may find T t such that for all t' >t>T e , 



< e, 



Jo 00 v(*)<fo(*)<k 



A" 



where 



Jo <t>d z ) dz 

K:=1 + (tP x , v)~ 2 



/ °V WaQ*)^ 
Jo°° M z ) dz 

^ s (P s p^lu)ds 



< e, 



Then, for t' >t>T e 



((99,^)-(^,X, 



= E,[((^,X t ) 2 ] +E,[(^,X^) 2 ] -2E,[(^,X t )( V 9,X i ,)] 
< 4e. 

Thus (ip, Xt) is a Cauchy sequence in L 2 and converges to a limit. In particu- 
lar, there is a random finite measure X^ such that Xt — > X^ in probability 
in the topology of weak convergence of finite measures on M + . 

It remains only to show that the limit is a multiple of the finite measure 
with density <j)\. For any bounded test function ip, let 



A t := (^X t )-(l,X t )(^^). 
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By (H3D, 



E[A 2 t ] 




(Pt<P,v) 
+ (if, cPx} 2 + {(p, 0a) 2 



Pt-sV 




p- 



ds, v 



Pt-sl ' 

(PM. 



ds, v 



which converges to as t — > oo. 
7.3. Optimal growth conditions. 



Pt_ a l Pt- a <p 

(P t l,u) (P t l,u) 



ds, v 



□ 



Theorem 7.5. Suppose b is a fixed drift function, (3 is a fixed non- 
increasing, non-constant net reproduction function, and (o~ m ) is a sequence 

of diffusion functions. Write (P^ m ^) t >o for the semigroup associated with 
a m , b, and (3, assuming complete reflection at 0. Assume that (|18p holds for 
each m. Denote the corresponding asymptotic growth rate by X m > P^. 
Suppose that either of the following two conditions hold. 



inf inf 

m a;GR+ 



b(x) 
o- m (x) 



and 



lim inf a m (x) 



=: 7 > 



oo, 



-oo 



(b) 



inf b(x) =: 6* > 



and 



Then lim^oo A m 



lim sup cr m (x) = 0. 



Po 



Proof. The condition and the conclusion remain true if (3 is changed by 
a constant, so we may assume, without loss of generality, that (3 is non- 
positive. This means that P t {m) is the semigroup of a killed diffusion Y^ m \ 
Let Tg be the killing time of Y^ m \ Write y( m ) for the diffusion without 
killing that corresponds to Y^ m \ Assume that Poo > —oo. The proofs for 
Poo = — oo are similar and are left to the reader. Since A m > Poo always 
holds, it suffices in both cases to establish that liminfm^oo A m < p^. 

Consider first the case of condition (a). It follows from the assump- 
tion that is a regular boundary point that the Liouville transformation 
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F m (x) := Jq du/a m (u) is finite for all x E M + . The transformed diffusion 
y(m) ._ p y-(m) ^ as diffogfon function the constant 1, drift function 



(24) 



a o F„ 



and killing rate function — /3 o F^ 1 . The asymptotic growth rate for 

is also A m . Write := F m o y( m ) for the diffusion without killing that 

corresponds to V^ m \ 

Since limm^oo inf 3 , gK + a m (x) = oo, it follows that linim^oo F~ 1 (x) = oo 
for all x > 0. Hence, given e > there exists M such that for m > M, 

/3of;V)<AW, 



where 

'0, 0<x<e, 
(/?oo + e)^, e<x<2e, 
k A» + e, a; > 2e. 

Write W for Brownian motion with drift 7 and complete reflection at 0. 
By the comparison theorem for one-dimensional diffusions and the assump- 
tion that (3 is non-increasing, 
r-t 



E 



cxp 



PoF-\vW)ds 







< E 



< E 



exp (J%oF- 1 (W s )ds 



cxp 



Pe(W s )d S 



if W m ) and W have the same initial distribution. Thus A m < A e , where 

A e := infjA : <fix,e( x ) i s non-negative for all x} 
for (fix e the unique solution to the initial value problem 



1 7 



(25) -</>'{, - 70' v + /3^A,e = A<AA, e with Ai£ (O) = 1 and A (O) = 7. 

Solutions of the equation 
1 



^f"(x)- 1 f'(x)-rf(x) = Xf(x) 

for some constant r are linear combinations of the two functions 
exp ( (7 - yV + 2(A + r)) x* 

and 



exp (^(7 + \/7 2 + 2(A + r; 
Solutions of the equation 



-f"(x) - 7 f(x) + {p- qx)f(x) = A/Or) 
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for some constants p and q are linear combinations of the two functions 

7 2 + 2(\-p + qxY 



e 7X Ai 
and 



2 2/3 g 2/3 
7 2 + 2(A - p + qx) 



2 2/3 g 2/3 

where Ai and Bi are the Airy functions (see Section 10.4 of [J). It follows 
that for fixed A the first and second derivatives of 4>\ ifi are bounded on [0, e] 
and of orders at most e _1//3 and e -2 / 3 , respectively, on [e,2e]. Thus 

(26) lim sup \4>\ ie (x) — 1| = and lim sup \(f>'\ e (x) — j\ = 0. 

40 0<2<2<e ' 40 0<x<2e 

Consider the initial value problem 

2$Ug - l<t>\,c,g + 0°c<t>\,c,g = ^4>\,c,g with 4>\, c ,g{°) = C aIld ^A,c, 3 ( ) = 9- 

Now 

inf{A : (f>x : i a (x) is non-negative for all x} 

is simply the negative of the asymptotic killing rate of reflected Brownian 
motion with drift 7 killed at constant rate —Poo, namely P^ (this also follows 
from explicit computation). If A > p^, then an explicit computation of the 
solution shows that 4>\, c ,g is non-negative for c sufficiently close to 1 and g 
sufficiently close to 7. 

Choose any A > /?oo- On (2e,oo), the eigenfunction coincides with 
one of the functions <fi\,c,g- From (|26p we know that by choosing e sufficiently 
small, we can make c as close as we like to 1 and g as close as we like to 7; 
hence, <j)\ te is non-negative on (2e, 00). The initial values of (|25p show that it 
is also non-negative on [0,2e] for e sufficiently small. Thus lim inf e |o Ae — 
since A was chosen arbitrary < P^, it follows liminf e | k e < Poo- Combining 
this with (jllh . we see that and the proof for condition (a) is complete. 

Now we consider the case of condition (b). Set 

Y? m = sup a 2 m {x) 



and 

z (m) = yim) _ Y (m) _ f* b( yM)ds _ ^(yM), 

Jo 

where £°(Y (m ^) is the local time at of Y^. Then z[ m) is a martingale 
with quadratic variation process bounded by S^t. Put 

4 m) * := sup \Z^\. 

0<s<t 

Because a continuous martingale is a time-change of Brownian motion, there 
is a universal constant C such that 

F{Z^>>e}<CeM-e 2 /2^ 2 m t). 
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For any 5 > 0, 
P{r^ >t}=E 

< E 



cxp 



cxp 



(3{Y^)ds 



St 



P(b*s - bj5/2)ds 



+ P{Z ( t m) * >b.t6} 



< exp I (1 - 5)t sup (3{x) ] + ¥{zl m> > KtS/2} 

V x>b„tS/2 J 

< exp((l - 5) p(b*tS/2) t) + C exp{-t(bJ) 2 /8T, 2 m ). 

(Recall we have assumed that (3 is non-increasing.) Since limm^oo S 2 
and lim^^oo (3{x) = (3^, we may find T such that 

P{r^ m) > *} < exp((l - 2<5)/3oot) 

for all t > T and all m sufficiently large. It follows that A m < (1 — 25)(3 O0 
for ?Ti sufficiently large. Thus lim inf m ^oo A m < [3^, as required. □ 







Appendix: Convergence of the branching processes 

We show in this section that the sequence of measure-valued processes 
(X n ) converges to a Dawson- Watanabe superprocess X under suitable con- 
ditions. We first recall the definition of the limit X as the solution of a 
martingale problem. 

For simplicity, we deal with the biologically most relevant case where there 
is no killing at in the damage and accumulation diffusion; that is, po = 1 
in the boundary condition (|8|), so that the domain D(£) of the operator 
L consists of twice differentiable functions <p that vanish at infinity, satisfy 
([8]) for po = 1, and are such that Lcf) is continuous and vanishes at infinity. 
Similar arguments handle the other cases. 

Suppose that the asymptotic rescaled branching rate p is bounded and 
continuous. Recall that the asymptotic net birth rate (3 is continuous and 
bounded above. Let M. denote the space of finite measures on M + equipped 
with the weak topology. For each v € A4 there is a unique-in-law, cadlag, 
.M-valued process X such that Xq = v and for every non- negative (b G D(£) 



exp(-{4>,X t )) - exp(-((f>,v}) 



exp(-(0,X s ))(-^ + ^ 



,X,) ds 



is a martingale. Moreover, X is a Markov process with continuous paths. 

Our convergence proof follows the proof of convergence of branching 
Markov processes to a Dawson- Watanabe superprocess, as found in Chap- 
ter 1 of [25] or in Chapter 9 of [28]. Because we only wish to indicate why 
a superprocess is a reasonable approximate model for a large population 
of fissioning organisms, we don't strive for maximal generality but instead 
impose assumptions that enable us to carry out the proof with a minimal 
amount of technical detail. 
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Let (ILJt>o be the Markovian semigroup generated by 

with boundary condition ([8]) for po = 1. That is, the generator of (ILjoo is 
the operator S acting on the domain D(S) consisting of twice differentiable 
functions (p that vanish at infinity, satisfy (JHJ) for pq = 1 with a 2 replaced 
by o"^ ot , and are such that S</> is continuous and vanishes at infinity. Thus 
(Ht)t>o is the Feller semigroup of the particle motion for the measure- valued 
processes X n . 

Let E denote the state space formed by taking the one-point compacti- 
fication of M + . The semigroup (Tlt)t>o can be extend to functions on E in 
a standard manner. The point at infinity is a trap that is never reached by 
particles starting elsewhere in E. The domain of the extended semigroup is 
the span of D(9) and the constants (which we will also call D(S)) and the 
generator is extended by taking 91 = 0. We also extend the measure- valued 
processes X n - the point at infinity never collects any mass if no mass starts 
there. 

Set I3 n {x) := a n {x) - k n {x) and p n {x) := n^a^x) + k n {x)). Suppose 
that j3 n and p n are bounded and continuous. Given a bounded continuous 
function g : R+ R, put 3+g(x) := / g(x + y)d-n™(y) and d-9(x) := 
f g(x — y) dTT™(y) 

We have for / G D(3) with inf^ f{x) > that 

exp((n log - exp((nlog/ ) X "}) 



exp((nlog/,Xf)) 
S/ + np n 



r + f-r.f + — -f 



np n np n 



,X2)ds 



is a cadlag martingale (cf. the discussion at the beginning of Section 9.4 of 
[28]). 

If Xq is a non-random measure v £ M, then essentially the same argu- 
ment as in the proof of Lemma 9.4.1 of [28] shows that 



(27) E„ [<l,Xf)] < (l,i/)exp \tsupp n (x)ir{x] 
and 

(28) p(sup(l,X")exp ( -t sup p n {x)/3 n {x) ) > zl < 



(the result in [28] has the analogue of sup^. \p n {x)(3 n {x)\ in place of 
sup x p n (x)P n (x), but an examination of the proof shows that it carries 
through with this constant). 

In order to motivate our next set of assumptions, imagine just for the mo- 
ment that the jump distribution 7r" is the distribution of a random variable 
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n 1 / 2 Y X where Y x has moments of all orders. For 4> sufficiently well-behaved 
and f n :=l — n~ l <f), a Taylor expansion shows that 
(29) 



fn(x) 

is of the form 



9f n (x) + np n (x) 



d n + fn(x)-d n .fn(x) + ^P,-f n (x) 



np n (x) np n (x) 



.2° 



.2 

mot 



(x)4>"(x) + b(x)(j)'(x) + I3 n {x)(j)(x) 



a n {x) + k n {x) 

plus lower order terms, and this converges pointwise to 

±a 2 (x)^"{x) + b(x)<f>'(x) + P(x)<Kx) - \ P {x)<t> 2 {x) 

under the assumptions (QJ) , ^ , and © . 

With this informal observation in mind, assume that 2)(£) C D(S) and 
that if f n := 1 — n~ l (f> with cj) G !>(£) and inf^ > 0, then the function 
in (f29|) converges uniformly to — L(f> + ^p^ 2 . 

Under this assumption, we can follow the proof of Theorem 9.4.3 in [28J 
(using (I28j) to verify the compact containment condition in the same man- 
ner that the similar bound (9.4.14) is used in [28]) to establish that if Xq 
converges in distribution as a random measure on E, then the process X n 
converges in distribution as a cadlag E- valued process to X. 
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